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By means of Monte Carlo techniques, we study the role of disorder on a system of hard-core 
bosons in a two-leg ladder with both intra-chain (t) and inter-chain (t') hoppings. We find that 
the phase diagram as a function of the boson density, disorder strength, and t! jt is far from being 
trivial. This contrasts the case of spin-less fermions where standard localization arguments apply and 
an Anderson-localized phase pervades the whole phase diagram. A compressible Bose-glass phase 
always intrudes between the Mott insulator with zero (or one) bosons per site and the superfluid 
that is stabilized for weak disorder. At half filling, there is a direct transition between a (gapped) 
rung-Mott insulator and a Bose glass, which is driven by exponentially rare regions where disorder 
is suppressed. Finally, by doping the rung-Mott insulator, a direct transition to the superfluid is 
possible only in the clean system, whereas the Mott phase is always surrounded by the a Bose glass 
when disorder is present. The phase diagram based on our numerical evidence is finally reported. 

PACS numbers: 05.30.Jp, 71.27.+a,71.30.+h 



I. INTRODUCTION 

Interacting bosons in one-dimensional (ID) or quasi- 
1D lattices are of interest in many physical contexts, 
ranging from Josephson-j unction arrays^ to more recent 
experiments on ultracold bosons loaded in ID optical 
trapsP Especially the latter ones offer the unique op- 
portunity to fine tune the experimental parameters and 
realize in laboratory a wide variety of bosonic models 
where kinetic energy, inter-particle interaction and dis- 
order can be varied at will. In particular, it is possi- 
ble to tune interaction to such an extend that atoms es- 
sentially behave as hard-core bosons confined along ID 
tubes witlPand without^ a superimposed optical lattice. 
In this condition, the huge repulsion prevents bosons 
from occupying the same position in space and induces 
a sort of Pauli exclusion principle. It is well known that 
hard-core bosons can be mapped onto spin-less fermions 
by the so-called Jordan- Wigner transformation. 5 How- 
ever, since the Jordan- Wigner fermions are non local in 
terms of the original bosonic operators, the local bosonic 
Hamiltonian is transformed into a very complicated non- 
local interacting fermion model. This approach simplifies 
in ID tight-binding models with only nearest-neighbor 
hopping: here, non-locality is absent and the bosonic 
model maps onto non-interacting fermions, easily solv- 
able. Whenever the lattice is not rigorously ID, i.e., when 
more chains are coupled together or longer-range hop- 
pings are considered, the corresponding fcrmionic prob- 
lem contains complicated interaction terms that may be- 
come highly non-local when the number of chains in- 
creases. 

The difference between hard-core bosons and spin-less 
fermions is even more pronounced in the presence of dis- 
order. Indeed, non-interacting spin-less fermions in a 



quasi- ID system are always Anderson localized for any 
disorder strength. On the contrary, hard-core bosons can 
become superfluid in the presence of disorder as soon as 
they can exchange among each other, a property that oc- 
curs already in the simplest case of a two-leg ladderP^ 
From this point of view, ladders of hard-core bosons rep- 
resent an ideal case study to uncover the role of Bose 
statistics versus Fermi statistics in the presence of disor- 
der. 

From the purely theoretical side, we think this is an in- 
teresting issue. Single-particle wave functions are always 
localized in a quasi-lD disordered lattice. It follows that 
any Slater determinant built with such wave functions is 
localized too, so any many-body wave function for non- 
interacting fermions. On the contrary, hard-core bosons, 
which like spin-less fermions cannot occupy the same site 
but whose wave function is symmetric under exchanging 
two particles, can cooperatively act and give rise to a 
delocalized superfluid phase. 

Also from the experimental side ladder systems ar e of 
interest, as they can be realized with optical lattices!^ 
In realistic experimental setups, a two-leg ladder can be 
realized through a double- well potential along a direction 
(say, y) like in Ref. [TUl and a potential creating a cigar 
geometry in the re-axis. Further superimposing a periodic 
potential along x, one could finally realize a two-leg Bosc- 
Hubbard model with tunable hopping rates among and 
between the legs. Disorder can be introduced by superim- 
posing a disordering lattice or introducing a speckle po- 
tential. We further mention that bosonic ladder systems 
are realized also in magnetic materials.^ For example, 
the disorder-free compound IPA-CUCI3 has been found 
to be a prototypical 5 = 1/2 antifcrromagnctic spin lad- 
der material, which can be thought as a system of inter- 
acting hard-core bosonsP^ Here disorder is introduced 
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by means of random chemical substitution, i.e., IPA- 
Cu(Cli_ l; Br a .)3. Neutron scattering experiments have 
shown convincing evidences of the spin-analogous of a 
Bose-glassP^ 

Given its both theoretical and experimental interest, 
we decided to study the phase diagram of a simple model 
of hard-core bosons hopping on a two-leg ladder with 
bounded on-site disorder by means of Green's function 
Monte CarloP^In short, we find that the phase diagram 
as function of the density and the ratio between inter- and 
intra-chain hopping includes three phases: a localized 
Bose-glass, a superfluid and, at half-filling, a so-called 
rung-Mott insulator that seems to be always surrounded 
by the glass. 

The paper is organized as follows: in Sec. [TTJ we intro- 
duce the model and the numerical methods; in Sec. |III| 
we briefly discuss the clean system; in Sec. |IV| we present 
our results for the disordered model; in Sec.|V| we finally 
draw our conclusions. 



II. THE MODEL 

We shall consider a system of disordered hard-core 
bosons on a L = 2 x L x lattice. The Hamiltonian reads 



i,rj— 1,2 i 
+ ^2 e i,V n i,V> (!) 



where b\ (b i ri ) creates (destroys) a boson at rung i on 
the chain r\ = 1,2. The matrix elements t and t' are the 
hopping amplitudes along legs and rungs, respectively. 
The disordered potential couples to the density operator 
n i.r) = b\ bi iV and it is described by random variables 
that are uniformly distributed in [—A, A]. Finally, the 
hard-core constraint of < 1 is implied. 

We study the Hamiltonian Eq. by Green's func- 
tion Monte Carlo with a fixed number M of bosons on L 
sites P^n = Al/L being the average density. The Green's 
function Monte Carlo approach is based on a stochas- 
tic implementation of the power method technique that 
allows, in principle, to extract the actual ground state 
I^Gs) °f a given Hamiltonian "H, from any starting trial 
(e.g. variational) wave function l^y), provided that 
(^vI^gs) 7^ 0. In order to improve the numerical ef- 
ficiency, it is important to consider a good starting wave 
function, for which we use one- and two-body Jastrow 
factors applied to a state where all bosons are condensed 
at momentum q = 0. The one-body Jastrow factor makes 
it possible to vary the local density of the bosons. On the 
contrary, the two-body Jastrow is taken to be translation- 
ally invariant.^ 








t'/t 








RMI 

/ 


3 










2 






SF 


SF 


1 


Metallic Phase 



0.25 0.5 0.75 n 0.25 0.5 0.75 n 

FIG. 1: (Color on-line) Upper panel: Superfluid stiffness as 
function of the density n at fixed t'/t — 2 on a lattice with 
L = 2 x 50. Lower panel: clean phase diagram of hard-core 
bosons (a) and spin-less fermions (b) on the two-leg ladder 
as function of t'/t and density n. The rung-Mott insulator is 
denoted by RMI, the superfluid by SF and the band insulator 
by BL 



III. THE CLEAN SYSTEM 

Before considering the disordered case, it is useful to 
briefly discuss the clean system, where = (see also 
Appendix). First, we consider the limit t' = 0, i.e., two 
uncoupled chains. In this situation, the ground state is 
a superfluid with quasi-long-range order for any density 
< n < 1. At densities n = and 1, there is a (triv- 
ial) "frozen" Mott insulator due to the infinite on-site re- 
pulsion, which completely suppresses charge fluctuations. 
Let us now analyze the the opposite limit t'/t^> 1. Ex- 
actly at half filling, i.e., n = 0.5, there is one boson per 
rung and the wave function can be approximately written 
as a independent product of single-particle rung states as 



|* GS )^JJ(6t 1+ &t 2 )|o). 



(2) 



The system is in the so-called rung-Mott insulator with 
a unique ground state and a gap to all excitations.^ At 
half filling, the transition between the rung-Mott insu- 
lator and the superfluid takes place exactly at t' = 0, 
since the inter-chain hopping represents a relevant per- 
turbation that immediately opens a gap in the excita- 
tion spectrum, see Appendix. This transition is of the 
Berezinskii-Kosterlitz-Thouless type, which makes it dif- 
ficult to observe by numerical simulations on finite clus- 
tersP 

Now, if few bosons are added or removed to the insu- 
lating state, a superfluid is stabilized. Therefore, at any 
other filling < n < 0.5 and 0.5 < n < 1 the system will 
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FIG. 2: (Color on-line) Upper left panel: superfluid stiffness 
p s as a function of the disorder strength A/t for different n 
and t! jt — 1. Upper right panel: superfluid stiffness p 3 as 
a function of the density for different A/t and fixed t' jt = 
2. Lower panel: low-density phase diagram of the hard-core 
bosonic model for t' jt = 1 and t! ft — 2. Calculations have 
been done on a L — 2 x 50 system. 




FIG. 3: (color on-line) Superfluid stiffness p s as a function 
of t'/t for different A/t on a two-leg ladder with L = 2 x 50 
sites. The density has been fixed to n = 0.4. 



be always superfluid, as can be seen in Fig. [T] where the 
superfluid stiffness as function of the density of particles 
n is shown for fixed t'/t = 2. The superfluid stiffness 
starts to grow at low density, it reaches a maximum, and 
finally vanishes at n = 0.5 at the rung-Mott insulator. 

We conclude this section by mentioning that, for spin- 
less fermions at half filling, there is a transition from a 
metallic to a band insulator at t'/t = 2, where a gap 
opens up. For all the other densities n ^ and 1, the 
ground state is metallic. The clean phase diagrams for 
hard-core bosons and spin-less fermions on a two-leg lad- 
der are sketched in Fig. [I] for comparison. 



IV. THE DISORDERED SYSTEM 
A. Low-density phase diagram 

Now we turn to the disordered case. Let us start by the 
low-density regime of the phase diagram at fixed inter- 
chain hopping as function of the density and disorder. In 
Fig. [2] we report our results of the superfluid stiffness as 
a function of disorder strength (for t'/t = 1) and density 
(for t'/t = 2); the low-density phase diagram is reported 
as well. We find that, for any finite disorder A/t, the low- 
density phase is a Bose glass that turns superfluid above 
a critical density. The trivial Mott insulator with zero (or 
one) bosons per site is therefore always separated from 
the superfluid by the Bose-glass phase. This is a remark- 
able result since, in a single chain with nearest-neighbor 
hopping only, hard-core bosons are equivalent to spin- 
less fermions, which Anderson localize for any density. 
Hence, in a two-leg ladder, hard-core bosons behave dif- 
ferently from spin- less fermions; while the latter ones re- 
main always localized, the former ones show a superfluid 
phase stabilized by the inter-chain hopping, as it was 
predicted using bosonization and renormalization-group 
techniques.^ The idea is that, in a strictly ID geometry 
with only nearest-neighbor hopping, the statistics of the 
particles does not matter. However, whenever particles 
may be interchanged (by non-strictly ID paths) bosons 
can form a superfluid, even in presence of disorder. We 
just mention that the same behavior holds also on a sin- 
gle chain with longer-range hopping. This scenario can 
be understood in very simple terms as follows. At very 
low fillings, the statistics of the particles does not matter 
so much and hard-core bosons, as free fermions, local- 
ize due to disorder, giving rise to the Bose glass. This 
is because the length over which the single-particle wave 
function extends is short enough that the wave functions 
of two particles never overlap. As soon as the filling is 
increased, the particles get closer to each other and the 
single-particle wave functions begin to overlap. At this 
point the statistics of the particles starts to play a role. 
If the particles are fermions they will still be localized 
(in D < 2), whereas bosons may stabilize a superfluid, 
as confirmed by our numerical simulations. 

In the range of values of inter-chain hopping that we 
have studied, the effect of a larger t'/t in the low-density 
phase diagram is to slightly reduce the superfluid re- 
sponse of the system, as can be also seen in Fig. [2] Al- 
though the actual thermodynamic value of the transition 
between the Bose glass and the superfluid may be rather 
different from the one obtained by our calculations, be- 
cause of strong size effects, the present results give a qual- 
itative correct insight into the phase diagram. Finally, 
we would like to mention that the exact behavior of the 
transition line between Bose glass and superfluid phases 
A c (n) cx n a is hard to be found by numerical calcula- 
tions. Although an almost linear fit is found, i.e., a = 1, 
a different power-law cannot be excluded, as implied by 
the arguments of Ref. [OH 
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FIG. 4: (Color on-line) Distribution P(E g ) of the gap as 
function of disorder strength and fixed density n — 0.5 and 
t'/t = 2 on a two-leg ladder with L = 2 x 50 sites. The clean 
gap is shown for comparison in the upper-left box as a blue 
bar. 

B. The effect of the inter-chain hopping 

Here, we want to investigate the effect of the coupling 
if on the otherwise insulating (Anderson localized) de- 
coupled chains. Generally speaking, for any value of dis- 
order, a certain finite ratio t'/t is necessary to drive the 
system into a superfluid phase, hence the system remains 
in the Bose-glass phase for small t' jt. However, for small 
disorder, the localization length of the Bose glass is ex- 
pected to be very large. This means that on clusters that 
are accessible to numerical simulations it may be very 
hard to see the Bose glass region. This fact is indeed 
confirmed by our results on the superfluid stiffness as a 
function of the inter-chain hopping, see Fig. [3} Rapidly, 
as a small t' /t is introduced, a large superfluid response 
is found for small disorder (e.g., A/t = 0.6 and 1.0 in 
the figure). It is also observed that, by a further in- 
crease of the inter-chain hopping the superfluid stiffness 
reaches a maximum and then eventually decays, since for 
t/t' the system decouples in a collection of decoupled 
rungs (which are obviously not superfluid). Moreover, as 
the disorder is increased, the superfluid response is sup- 
pressed, until the system cannot attain superfluidity any 
longer and remains localized for any value of t' . 

C. The rung-Mott phase in presence of disorder 

The effect of disorder on the rung-Mott phase at den- 
sity n — 0.5 system is now discussed. From general 
grounds, it is expected that the presence of disorder will 
fill the gap with localized states, so to induce a transition 
to a gapless Bose-glass phase. In practice, given a ratio 
t! /t, the rung-Mott insulator will survive up to a certain 
critical value of A/t, where the gap will be completely 
filled and the system will become compressible. This sit- 
uation is similar to the one of the Bose-Hubbard model 
at integer fillings^, where a direct transition between the 
Mott insulator and the Bose-glass phase is expected by 
decreasing the ratio between the on-site repulsion U and 



the disorder strength ApStHl Also in our case of hard- 
core bosons, we can make use of the argument based on 
the fact that, if A is larger than half of the energy gap 
of the clean insulator Eg lean , then the ground state must 
be compressible; otherwise the system is incompressible 
with a reduced gap given by E g — E c g lcan — 2 A. In partic- 
ular, for t'/t > 1 we have that Ef ean ~ 2t' . Therefore, 
the gap will vanish around a critical value of disorder 
A c ~ t' . These arguments should hold exactly only in 
the infinite system and large size effects are expected be- 
cause this transition is of the Griffiths type, i.e., driven 
by exponentially rare regions which are locally ordered.^ 
On the other hand, on any finite system the transition 
from the gapped to the compressible phase will appear at 
a larger A c , since these exponentially rare configurations 
will be hardly sampled on finite clusters. 

We recently proposecP^ a method to alleviate the 
strong size effects that consists of computing directly the 
distribution probability of the gap 

P(E g )=Y / t(E g -^+^), (3) 

where fi^ — ± (£^ ±1 — Efo) (E^ being the ground- 
state energy with M particles on the realization a of 
disorder). This definition of the gap distribution is in- 
troduced because in disordered systems the gap can be 
overcome by transferring particles between two rare re- 
gions with almost flat disorder shifting the local chemical 
potential upward and downward. These exponentially 
rare regions may be far apart in space and represent rare 
fluctuations (Lifshitz's tail regions), thus it is useful to 
imagine that a large system is made by several subsys- 
tems, each represented by a different disorder realization 
of our L-site cluster, and construct the gap by using the 
process of taking one particle from region a to region 
/3p2uf such processes are allowed at no energy cost, i.e., 
P(0) 0, the corresponding system will be gapless. One 
could define an alternative estimate of the gap as 

E™ = mm at0 \»+- fJ ,- p \, (4) 

with all the disorder realizations a and /3. 

In Fig. |1J we show the distribution probability at half 
filling as function of disorder, t'/t = 2 and L = 2 x 50. 
For small values of disorder the gap survives, while for 
A/t = 2 the probability to find zero gap is finite, which 
we interpret as signalling zero gap in the infinite system 
and a Bose glass phase. We note that, when considering 
the case of the rung-Mott insulator, this method per- 
forms a bit worse than in the case of the Bose-Hubbard 
model (at integer fillings). Indeed, for this value of the 
hopping parameters, we have that Eg lean /t ~ 1.11, giv- 
ing rise to A c /t ~ 0.55, which is much smaller than the 
value obtained by numerical simulations. However, we 
would like to mention that, even though the finite-size 
analysis of P(E g ) overestimates the actual value of the 
transition, it gives a sizable improvement with respect to 
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FIG. 5: (Color on-line) Left panel: Superfluid stiffness p s as 
a function of the density and several values of A/t and fixed 
t'/t — 2. The stiffness of the clean system is also shown for 
comparison. Right panel: Superfluid stiffness p s as function 
of the disorder bound for several densities close to half filling 
and fixed L = 2 x 50 and t'/t = 10. 




FIG. 6: (Color on-line) Results of the 2 x 50 lattice for 
the phase diagram in the vicinity of the rung-Mott phase for 
t'/t = 10. 



the simple calculation of the average gap value in pres- 
ence of disorder. The precise determination of the criti- 
cal point is well beyond any numerical calculations and, 
therefore, we take advantage of the fact that its estima- 
tion can be done by using the criterion of Ref. 
namely A c = Ef can /2. 



D. Transition from the superfluid to the gapped 
phase 



to the rung-Mott insulator as the density n — ¥ 0.5, or 
there is always an intruding Bose glass phase. In this 
regard, whenever the gapped state is doped with a few 
particles or holes such that their typical spacing will be 
large, those few carriers on top of the rung-Mott phase 
will effectively see a disordered background. Therefore, 
standard single-particle Anderson localization arguments 
apply and the system remains insulating by localizing 
those few carriers on the Lifshitz's tails that are in the 
Mott gap. By further increasing the density of particles 
(or holes), a superfluid is eventually formed. This sim- 
ple single-particle argument implies the presence of an 
intervening Bose glass between the rung-Mott phase and 
the superfluid. We proceed to test this argument quan- 
titatively. In Fig. [5j we report our numerical results for 
the superfluid stiffness at densities close to n = 0.5 and 
t'/t = 2. Our data is consistent with a transition driven 
by density from the superfluid phase through the Bose 
glass to finally end up with the rung-Mott insulator. For 
example, for A/t = 1, the superfluid stiffness appears to 
vanish just before n = 0.5, however the region in which 
the Bose glass takes place is very small. For A/t = 2 the 
rung-Mott insulator has already been wiped out by the 
effect of disorder, as observed in Fig. [4j and, therefore, 
this issue cannot be addressed. 

By considering a larger value of t'/t we have two ad- 
vantages: first, the Mott gap is larger, such that the 
gapped phase is more robust, second, the localization 
due to disorder is expected to be enhanced, thus, enlarg- 
ing the Bose glass region at n < 0.5 (or n > 0.5). These 
facts enable us to provide further evidence in favor of 
an intervening Bose glass in between the rung-Mott and 
the superfluid as follows. We have performed simulations 
with a rather large t'/t — 10. For such a value of the hop- 
pings, we have that £!| lean ~ 17.31, such that the tran- 
sition from the gapped to the compressible (Bose-glass) 
phase is argued to occur at a A c ~ 8.65. Therefore, the 
system is expected to be gapped for A < 8.65. In Fig. [5] 
we present our results for the superfluid stiffness as func- 
tion of the disorder strength for several densities close 
to half filling. From these calculations we can easily see 
that, as the density approaches n = 0.5, the critical point 
where the stiffness vanishes gets smaller, leaving room for 
a large Bose-glass phase in between the superfluid and 
the rung-Mott insulator. Given our results, we can draw 
the phase diagram for densities close to the rung-Mott 
phase, see Fig. 6] Notice that the large value of t'/t (and 
therefore the clean Mott gap) ensures the existence of a 
truly gapped state at half filling. All together, we can 
make the safe statement that the transition between the 
superfluid and the rung-Mott phases is not direct, but 
through an intervening Bose-glass state. 



In connection with the rung-Mott phase at half filling, 
we investigate the phase diagram in the vicinity of such 
a phase and address the question of whether it is possi- 
ble to have a direct transition from the superfluid phase 



V. CONCLUSIONS 

We have studied hard-core bosons on disordered two- 
leg ladders by using both numerical techniques and an- 
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FIG. 7: (Color on-line) a) Tentative zero-temperature phase 
diagram of hard-core bosons on the two-leg ladder system, b) 
Same phase diagram in a) but for a larger value of disorder 
A/t. 



alytical arguments borrowed from similar problems in 
Bose-Hubbard models. We have shown that the zero- 
temperature phase diagram is rather rich and contains 
different phases; apart from the trivial Mott insulators 
at n — and 1, that are totally frozen due to the hard- 
core constraint, we found superfluid, Bose-glass, and 
rung-Mott phases. This contrasts the case of spin-less 
fermions, where no metallic phases are possible and An- 
derson localization takes place for any density n =/= 0.5 
at finite A. A final sketched phase diagram, based upon 
our results, is reported in Fig. [7] In the case of no disor- 
der, i.e., A = 0, the superfluid phase pervades the phase 
diagram for all densities < n < 0.5 and 0.5 < n < 1 
and all t'/t ^ oo. When considering a finite disorder 
strength, the superfluid shrinks and a Bose-glass phase 
appears. Most importantly, the transition between the 
Mott and the rung-Mott phases and the superfluid ones 
is never direct, like in the Bose-Hubbard model. 20 23 
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Appendix A: Basic bosonization formulas 

There are several papers that discuss at length the 
harmonic-fluid representation of bosonic lattice Hamil- 
tonians following the seminal work by HaldanePf Nev- 
ertheless, we believe that it is worth listing some useful 
formulas, referring the interested readers to existing lit- 



erature for further detailsf^H^Zl 

In the long- wavelength limit, the boson density and 
creation operator on each chain a — 1,2 can be written 

as 



p a (x) 



Po 



i2ml tp a (x)-\-7rpQX 



m—~oo 



OO 

E 



(x) 



where po is the average density and the two fields 4>{x) 
and 6{x) satisfy 



4> a {x),S79p{y) = iir5 a/3 5(x-y). 



(Al) 



In the case of hard-core bosons, it could be useful to 
define (j>{x) and 6{x) in terms of right (R) and left (L) 
chiral fieldsP 

<f>a{x) = ^(j> a Il(x) +4*aL{x)^j, (A2) 

= * (<Wz) - <W^)) • (A3) 

The dimension d of each operator A(x), defined through 
(A(a;) A(0)) ~ x~ 2d , can be easily evaluated by recalling 
that 



/ e i fl>R( x ) e -*70n(O)\ 



l e il<pL{x) e -ij(j>L{0)\ 



In the two-leg ladder it is convenient to introduce the 
symmetric (s) and anti-symmetric (a) combinations <j) s = 

(<t>i +02^ /V% and (j) a — — /V%, respectively (and 
seemingly for 9 S and 9 a ). It follows that the inter-chain 
hopping becomes 



ip\(x)ip 2 (x) 



e . w -~ay~, e i2irpa(m+n)x x 



x gi\/5 (m+n) <f>, (x) giV? (m—n) <j> a (x) 



The first term in the right hand side is the most rele- 
vant one and opens a gap in the anti-symmetric sector 
such that 9 a acquires a finite average value while 4> a has 
exponentially decaying correlations. It follows that the 
leading operator generated by the inter-chain hopping is 



ip\{x) ^ 2 0*0 



-i\/28 a (x) 



1+2 cos 



: (x)+-llTpoX 



. (A4) 

In other words, t' not only gaps the anti-symmetric sec- 
tor but also generates an umklapp scattering in the sym- 
metric channel that becomes marginally relevant at half 
filling, where 4npo = 2tt. It is just this umklapp that is 
responsible for the appearance of the rung-Mott insulator 
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at any finite t' <C t. In addition, being marginally rele- 
vant, it opens a gap in a Berezinskii-Kosterlitz-Thouless 
fashion, which is hard to detect numerically. Following 
Ref. one finds that disorder gives rise to a fluctuating 
umklapp 



"Humklapp = / dx£(x) COS ^8 (j) s (x)j . 



(A5) 



with £,(x)£(y) = u 2 5(x — y) in case of a Gaussian noise. 
Therefore, at half filling the coupling constant of the 
umklapp has a finite value plus a fluctuating one £,(x); 
therefore, if the latter one is small, the gap is on aver- 
age finite (the Mott phase), while, for u above a certain 
threshold, the gap is washed out by disorder, leading to 
the Bose glass. 

Away from half filling, we should keep into account 
a renormalization of the symmetric sector that can 

be par ameteriz ed by a Luttinger liquid parameter K s 
througrPMEU 



Since the full density 
p(x) = p 1 (x)+p 2 (x) 



■V&,(a;) 



V2K; 



where L x 
Po ¥= 1/2 



V0 s (x), 



is the number of sites per chain. When 
the umklapp scattering in Eq. ( |A4[ ) ceases 
to play a role and what survives is just the disorder- 
generated term (A5). Conventional scaling arguments 
predict that such a term is relevant when K s < 3/4p3 
This would suggest that for K s > 3/4 a superfluid phase 
is stable, otherwise disorder is relevant and the Bose glass 
occurs. According to the theory of the commensurate- 
incommensurate transition in ID, we expect in the clean 
case that, as the density po approaches half filling K s — > 
1/2, which would imply that disorder becomes relevant 
already before the Mott transition is approached in den- 
sity. However, even if the density po — > we should 
expect K s — » 1/2. Therefore, whatever is the behavior 
close to half filling, it must be qualitatively the same also 
close to zero filling. We know that, at very low density, 
bosons localize in the Lifshitz's tails and superfluidity 
can arise only when the localization length becomes of 
the order of the interparticle distance. This also suggests 
that a harmonic-fluid representation is likely inadequate 
at low density, hence that the scaling criterium K s > 3/4 
for the appearance of superfluidity may not work. Seem- 
ingly, the same argument must apply close to half filling, 
so that it must not be surprising that the phase boundary 
between the superfluid and the Bose glass goes smoothly 
and almost linearly to zero, see Fig. [5j 



K s can be easily extracted by the static structure factor 
in momentum space 



(p(q) p(-q)) = 2K t 



qL x 

2tt ' 
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